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Abstract 

Evolution of the parton densities at NLO in as using = Q'^{1 — z)/z instead 
of the usual for the scale of the running coupling as is investigated. While this 
renormalisation scale change was originally proposed as the relevant one for x — > 1, we 
explore the consequences for all x with this choice. While it leads to no improvement to 
the description of DIS data, the nature of the gluon at low x, low is different, avoiding 
the need for a 'valence-like' gluon. 



1. Motivation — Non-singlet evolution 

In deep inelastic scattering (DIS) the dependence of flavour non-singlet quantities is 
quite straightforward. Taking a non-singlet structure function, e.g. F^g = — Fg^" or 



NS 



then at next-to-leading order (NLO) we have 
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where t = InQ^ , q{x,t) = xq{x,t), the relevant combination of quark and antiquark densi- 
ties (weighted by the appropriate charges squared) is denoted by g^vs and B^{z) is the quark 
coefficient function given, for example, in the MS scheme by 
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is the well-known LO q-q splitting function. 

Here [...]+ denotes the standard regularised functions defined by 

'dz [f{z)]^g{z)= f'dzfiz) [9{z)-g{l)] 
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The NLO evolution of the non-singlet quark density is governed by the qq and qq splitting 
functions 
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The NLO splitting functions take the form 

= ClPF{z) + h0FCAPG{z) + ]^CFNFPMAz) 

pS^{z) = iCl-lcFCA)PA{z) (6) 

and the explicit expressions for Pf{z), Pg{z), Pn^^z) and Pa{,z) can be found in for example. 

By combining eqs.(|I],^), the evolution of the non-singlet structure functions to 0{a'g) may 
then be expressed in the form 
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If, in the above as is a function of t only, the running couphng can be taken outside the 
convolution integral and we have the usual NLO in ^5 evolution as implemented in standard 
analyses of DIS. 

Next consider the case where the relevant scale depends on z as well i.e. asiQ'^) 
as{4>{z)Q'^) and, in particular, the choice = {1 — z)/z. The quantity W'^ = (5^(1 — z)/z is 
just the com energy squared of the virtual photon-parton scattering which controls the maxi- 
mum tranverse momentum occuring in the ladder graphs which are summed to give the leading 
log contribution. As 2; — >• 1, W"^ as well as is large and it has been argued p, ^ that, in this 
region, the relevant variable to account for the large logs which arise beyond the control of the 
renormalisation group is W"^ or a quantity closely related. 

Expanding as(t + ln[(f){z)]) to 0(q;|) we get 

^^^^^^^ = ^-Wn[4>U)](^J (8) 

which means that, to this order, the change of scale of as is equivalent to adding an extra term 
inside the [...]+ of eq.(^ equal to 

-kPo(^^Jp^^\z)\nm] (9) 

The above expressions for the evolution for and F^g are now to be understood with 
as = asit,z) but we can see from eqs.(|^,H) that to the shift in scale is equivalent to 

a term in Bg{z) equal to 2P^g\z)ln[(f){z)] or a term in P^l\z) equal to -|/3oP,(°)(z) ln[0(z)]. 
Thus taking (j){z) = {1 — z)/z the first term in eq.(|^) is generated by the change of scale 
Q"^ ^ Q'^{1 — z)/z and then to avoid double counting at we should use the simpler 

form 

B,{z)^ -|pW(^) + i(9 + 5z) (10) 

Since Po = ^Cf{11Ca — '^Np) we note that some of the Np dependence in the combination 
of the NLO NS splitting function and the quark coefficient function has been absorbed by the 
change of scale. One can go further and compute the form (j){z) should take in order to generate 
the entire Np dependence at Wong0 showed that can be achieved by taking 

0W - 0(^) = (^) exp[{Cp/P^^\z))C-^) - f|] (11) 

This is the BLM procedure Q where the Np dependence in a particular process is identified 
as arising from the vacuum polarisation contributions to that process and these are then entirely 
absorbed into the running coupling thus providing a method of summing such contributions 
to all orders. This procedure is therefore quite attractive but it requires the choice of the z- 
dependent scale to vary from process to process. In that sense the evolution of the non-singlet 
structure function is a different process from the singlet one. It is impractical to attempt 
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an analysis of DIS where the choice of renormahsation scale differs for the singlet and non- 
singlet combinations and here we shall explore the the consequences for a single choice of scale, 
(j){z) = {1 — z)/z. It is clear that this simple choice accounts, at (9(q;|), for large terms in the 
MS coefficient functions ocurring both in singlet and non-singlet expressions and it is worth 
investigating the phenomenological consequences of the potentially large logarithms which are 
summed by this procedure over the entire range of x. 

2. Quark singlet and gluon evolution 

Consider now the evolution of a singlet structure function, Fs{x,t) where 

Fs{x, t) = qs{x, t) + (^g) Bg{z) ® qs{x/z, t) + (^g) Bg{z) ® g{x/z, t) (12) 

where qs{x,t) = xJ2{(l{Xyt) + Qi^^t)), and Bg{z) is the gluon coefficient function given in the 
MS scheme by 

B,{z) = 2Pf^{z) In f + 2N,[%z{l - z) - I] (13) 



Sg 
where 

Pf^{z) = N,[z' + {l-zf] (14) 

is the LO q-g splitting function. 

The NLO evolution of the singlet quark density is governed by the appropriate SS and Sg 
splitting functions 
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where ^Pjjl^\z) and Ps'q{z) have the form 



APW(^) = CpNfF,,{z) 
P^sk^) = lCAN^Fil\z) + '^CfNfF^'\z) (16) 



and the relevant expressions for the F's can be read off from ref . . As in the non-singlet case, 
we can combine eqs.( P^JT3| ) and the evolution of the gluon into a 'one-step' evolution of the 
singlet structure function^ correct to 0(a;|) which is of the form 



(17) 

where Vg'\ V'g \ P'-"'^^ and S''^^ are given in ref 



Using eq.(||) we have again that the effect of changing the scale Q"^ Q"^ {1 — z) j z is to 
generate, at the logarithm terms which explicitly appear in both the MS gluon and 

quark coefficient functions. So in addition to eq.(|TO|) we take 

Bg{z) — > 2Nf[Sz{1 -z)-l] (18) 



For the gluon evolution we simply insert counter terms involving (j)[z) to restore eq.(|T7D at 
0{a\{t)) generated by the scale change. 
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While in the non-singlet case a single change of scale for as (eq.(|TT])) can be chosen to absorb 
all the Nf dependence at (9(a|), the situation is more complicated here and a similar BLM 
procedure would involve different choices of scale for the quark singlet and gluon evolution. An 
alternative to eq. (p!7D would be a one step evolution expressed in terms of physical structure 
functions only, Fs{x,t) and Fi^x^t) in which case the modifications due to — >• (l){z)Q'^ could 
be entirely absorbed at C(q;|) into the coefficient functions Bq{z) and Bl{z). 

3. Fitting the DIS data 

A practical problem that arises in using (f){z) = {1 — z)/z is simply that the integrations 
involved in the convolution integrals run up to z = 1 and so we must adopt some sensible 
approach for computing the running coupling at low values of the scale. Since the argument 
of as is actually timelike, then beyond leading order the coupling is complex. (Recall that 
the perturbation expansion is derived strictly in spacelike region and then continued analytically 
to the timelike region.) We consider two possibilities (i) compute asdVT^I) in terms of A^ in 
the standard NLO way but 'freeze' its value for < Q'j and (ii) use |as(W^^)| as the running 
coupling. This latter variable is claimed]^ to be a far more efficient expansion parameter (since 
it resums large terms involving vr^ arising from the analytic continuation) and grows only weakly 
as the argument becomes very small. The imaginary part of as is computed from the modulus 
at NLO, 



Im(a5) 



Poi 1 2 , /^l 

4 lovr 
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(20) 



Fig. displays the running couplings used in our fits compared with the 0:5 used by a 
conventional fit using a z-independent scale, the values of Aj^ giving the best fit to the DIS 
data being shown. For the choice asd^T^I) the quality of the fit is not overly sensitive to the 
precise value of the scale at which the coupling is frozen and we take a value around 1 GeV^. 
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As with the usual analysis of DIS structure functions, the parton distribution functions 
(pdf's) are parametrised at some starting scale, Qq, but now evolved according to eqs.(0,|15|,|T9|). 
In the recent MRST fits|^], the starting scale is Qq = 1 GeV^ and the gluon at small x is 
suppressed (so-called 'valence-like' gluon) at this in order to describe the relatively small 
slope dF2/dlnQ'^ of the structure function observed at HERAfl 



At 



lT| at low X and low 

low X, the mean value of z in the evolution is also low implying that the scale is much larger 
than Q"^ and hence the gluon evolution is naturally suppressed at low x. In Fig. ^ we show a 
comparison of the structure function evolution for a common set of pdf's at according to 
the two scales. The starting scale is Qq = 0.5 GeV^ and it is clear that evolving with the scale 
leads to dramatically slower evolution at low x and low Q^, even with a common value of 

The concept of a gluon distribution vanishing as x ^ (even at low Q"^) seems unnatural 
and leads to problems when attempting to evolve down in since physical quantities, such 
as Fl quickly become negative. The comparison in Fig. Q shows we can instead start with a 
larger gluon and in fact we find the data can be fitted a gluon distribution which is actually 
singular for Ql as low as 0.5 GeV^. 

DIS data fitted include the HERA data|[TD|, [IT], |g,BCDMS|[TB|, E6650, SLAC[|g] and 
NMC[ |T^] . We label the two types of fits as type (i) or type (ii) depending on the choice asdlV^I) 
or |a5(W^^)|. There has been much interest in the values of the slope dF2/d\nQ'^ shown by the 
low X, low 1995 ZEUS data[T^ that has prompted speculation about a possible breakdown 
of the standard theory |T^. A type (i) fit gives a satisfactory description of the slopes observed 



by ZEUS without the need to invoke a 'valence-like' gluon distribution at low — as shown 
in Fig. ^ and as anticipated above. Because of the more conventional behaviour of the gluon 
at small x one can evolve to values of below 1 GeV^. At the starting scale Qq = 0.5 GeV^ 
the gluon has the form 



xg{x, Q 



1.25x 



-0.046/ 



X 



,5.52, 



1 + 0.032v^ + 5.66a;) 



(21) 



which is virtually 'flat' at small x. 

Fig. ^ clearly shows a good description of the F2 slopes in the HERA range both for the 
conventional type of fit and for the new type (i) fit using a 2;-dependent scale. Copmparing in 
detail the quality of the fits, we obtain 





HI 


ZEUS 


BCDMS 


NMC 


SLAG 


E665 


No. data pts 


221 


216 


174 


130 


70 


53 


(MRST) 


164 


264 


249 


143 


116 


56 


(type (i) fit) 


164 


270 


243 


170 


134 


52 


(type (ii) fit) 


178 


298 


260 


223 


156 


48 



In the new fits, the cuts applied are the same as for MRST, in particular > 2 GeV^ 
except for the HERA data where we take > 1.5 GeV^. For the HERA data with Q"^ > 1.5 
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GeV^, fit (i) achieves virtually the same quality as the conventional type fit but can, in addition, 
give a good description down to Q"^ = 1 GeV^. The type (ii) is slightly worse for the HERA 
data. For the 'intermediate' x region covered by NMC, the new fits are not as successful as 
the conventional MRST fit, especially the type (ii) fit, and this is due to problem of trying to 
describe the slopes of the NMC data while simultaneously being consistent with the slopes 
measured at HERA. The MRST fit already underestimates the slopes observed by NMC and 
the new fits only serve to emphasise the disagreement as shown in Fig. ^. Type (ii) fit being 
even worse indicates that the results are sensitive to the precise nature of the coupling at very 
low values of the scale. 

4. Conclusions 

Although the motivation for evolving partons with a scale (f){z)Q'^ where (f){z) = {1 — z)/z 
stems from a procedure for resumming terms which are potentially large as x — >■ 1, it is worth 
exploring the phenomenological consequences for quarks and gluons over the whole x region. 
In fact we did examine whether the large x region was better described in terms of a running 
coupling with the choice of scale Q^(l — z)/z - the hope being that the observed strong 
observed at values of x beyond 0.6 might be absorbed by such a choice instead of conventional 



higher twist contribution ( see ref. [|19| for example) but we found no evidence for this. The 
change of scale (f){z)Q'^ generates corrections at C(a|) which already appear explicitly 

as the logarithm terms in the MS coefficient functions, so evolving with the 2;-dependent scale 
implies somewhat simpler expressions for the coefficient functions. It should be remembered 
that our change of scale is a change of renormalisation scheme and in order to maintain the 
expressions for the evolution of the physical F2 structure function to 0{al) we must modify the 
coefficient functions and the relevant splitting functions as shown in sections 1 and 2. Also the 
regularising counter terms (proportional to 5(1 — z)) are more complicated due to depending 
on z. 

Performing fits with two alternative prescriptions for handling the running coupling at low 
values of the scale leads to a preference for the simpler choice of freezing its value at around 1 
GeV^. In this case one can get acceptable fits to the DIS data but no improvement over the 
standard procedure is observed; in fact the problem of trying to reconcile the slopes dF2/d\nQ'^ 
measured by NMC with those measured at HERA is is actually aggravated. This subtle 
dependence of the slopes has only been successfully described when effects beyond standard 
DGLAP physics are included. In particular Thorne|]18[ was able to achieve a high degree 



of consistency between all the DIS data sets with a leading order, renormalisation scheme 
consistent calculation which included leading ln(l/x) terms. 

An advantage which has been gained in using the scale (5^(1 — z)/z is that the gluon 
distribution at low seems more 'natural' than the 'valence-like' gluon which the standard 
DGLAP description finds necessary to account for the small values of the slopes observed by 
ZEUS|jl2|. We therefore expect that quantities dominantly governed by the gluon are sensitive 



to the change of scale, in particular the longitudinal structure function at low x, low is quite 
different. Fig. |] shows a comparison between predictions for Fl at = 1 GeV^ resulting from 
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fits to F2 using or W"^ as the choice of scale. The curious 'dip' of the MRST curve is due 
to the suppressed gluon contribution at low x while the smoother behaviour of the fit with the 
z-dependent scale reflects a dominance of the gluon contribution at small x. 

In summary, the effects of resumming some of the log(l — x) terms through the change 
of scale does not lead to any improvement phenomenologically though it is consistent with a 
significantly different, and perhaps rather more natural, form of the gluon distribution at very 
low values of Q^. The exercise indicates that there is therefore a degree of uncertainty in the 
detailed nature of the gluon density at such low scales. 
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Figure 1: Running couplings used in the various fits. The sohd hue shows as'diy^l) versus W^, 
the dashed hne shows |q;s(H^^)| versus W'^, the values of Aj^ being indicated. The dotted hne 
shows, for the sake of comparison versus as{Q^) versus Q^, the running coupling used in the 
standard MRST fit. 
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Figure 2: Comparison of the evolutions with scales and from a common set of parton 
distributions at Qq = 0.5 GeV^ and with a common value of A.^^. 
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Figure 3: Slope of F2 for different x values measured by ZEUS 1 12]. Note the strong correlation 
between the value of x and the mean value of Q^. The data are compared to fits to DIS data 
using a scale or a scale W^. 
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Figure 4: Slope of F2 for different x values measured by NMC||TB||. Again the values of x and 
the mean value of are correlated. The data are compared to fits to DIS data using a scale 
or a scale W"^. 
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Figure 5: A comparison of Fl^x^Q"^ = 1.25) predicted by two fits, one (MRST) using the 
standard evolution scale, the other the type (i) fit using a scale W^. Also shown are the 
respective gluon distributions at the same — 1.25 GeV^ which are responsible for driving 
Fl. 
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